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(57) Abstract: Disclosed is a method for automatically establishing a system of equations according to the finite element process. 
Said system of equations relates to nodes of a structure comprising a body and an adjacent layer (Kl). Finite elements are created 
for the body. The layer (Kl) is automatically subdivided into volume elements by using geometrical data about the layer as well as 
specifications for meshing the layer. A finite element (1 10.2) of the body and one closest point (205.1) on said finite element are 
determined for at least one node (201.1) that is pait of a volume element of the layer, said finite element (1 10.2) being located closest 
to said node. A function is created for a physical relationship between the value that the physical variable takes at the closest point 
(205.1) and the values which said variable takes at the nodes of the closest finite element (100.2). The value of the physical variable 
at the determined point is eliminated by using the function when the system of equations is established. The invention allows the 
number of unknowns to be reduced in the system of equations. 



(57) Zusanunenfassung: Die Erfindung betrifft ein Verfahren zum automatischen Aufstellen eines Gleichungssystems gemUss der 
Finile-Elemenle-Methode. Das Gleichungssystcm beziehl sich auf Knotcnpunkte einer Konstruktion, die cinen K(3rper und eine 
W angrenzende Schicht (Kl) umfasst. Finite Elemente fUr den K5rper werden erzeugt. Die Schicht (Kl) wird automatisch 

1^ [Fortsetzung auf der ndchsten SeiteJ 



wo 2004/109559 A2 lillliiliilliiililllillili^ 



KG, KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, 
MG, MK, MN, MW, MX, MZ, NA, Nl, NO, NZ, OM, PG, 
PH, PL, PT, RO, RU, SC, SD, SE, SG, SK, SL, SY, TJ, TM, 
m TR, TT, TZ, UA, UG, US, UZ, VC, VN, YU, ZA, ZM, 

zw. 

(84) Bestiininungsstaaten (soweit nicht anders angegeben, fur 
jede verfUgbare regionale Schutzrechtsart)i ARIPO (BW, 
. GH, GM, KE, LS, MW, MZ, NA, SD, SL, SZ, TZ, UG, 
ZM, ZW), eurasisches (AM, AZ. BY, KG, KZ, MD, RU, 
TJ, TM), europaisches (AT, BE, BG, CH, CY, CZ, DE, DK, 
EE. ES, Fl, FR, GB, GR. HU, IE, IT, LU, MC, NL, PL. PT, 



RO, SE, SI, SK, TR), OAPI (BF, BJ, CF, CG. CI. CM. GA, 
GN, GQ. GW, ML, MR. NE, SN, TD. TG). 

VerofTentUcht: 

— ohne intemationalen Recherchenbericht und emeut zu ver- 
offentlichen nach Erhalt des Berichis 

Zur Erkldrung der Zweibuchstaben- Codes und der anderen Ab- 
kUrzungen wird auf die Erkldrungen ("Guidance Notes on Co- 
des and Abbreviations*') am Anfang jeder reguldren Ausgabe der 
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in Volumenelemente zerlegt, wobei geometrischen Informationen liber die Schicht und Vorgaben fur die Vemetzung der Schicht 
verwendet werden. Filr mindeslens einen Knotenpunkt (20L1), der zu einem Volumenelemeni der Schicht gehflrt, werden ein dem 
Knotenpunkt nMchstliegendes Finites Element (100.2) des KGrpers und ein nMchstliegender Punkt (205.1) auf diesem Finiten Element 
ennittelt. EineFunklion fUreinen physikaJischen Zusammenhang zwischendem Wert, den die physikalische Grtisse im nachstliegen- 
den Punkt (205.1) annimmt. und den Werten, den diese Gr&sse in den Knotenpunkten des nMchstliegenden Finiten Elements (100.2) 
annimmt, wird erzeugt. Beim Aufstellen des Gleichungssyslems wird der Wert der physikalische n Grdsse im ermittelten Punkt durch 
Binsetzen der Funktion climinicrt. Die Erfindung spart Unbekannte im Gleichungssyslem cin. 
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Verfahren zur Simulation einer Fugeverbindimg 

Die Erfindung betrifft ein Verfahren zum automat ischen Auf- 
stellen eines Gleichungssystems zur Beechreibung eines physi 
kalischen Verhaltens eines vorgegebenen Systems gemalS der Fi 
nite-Eleraente-Methode . 

Die Methode der Finiten Elemente ist aus ,,Dubbel - Taschen- 
buch fiir den Maschinenbau" , 20. Auflage, Springer- Verlag, 
2001, C 48 bis C 50, sowie aus T, R. Chandrupalta & A. D. Be 
legundu: „ Introduction to Finite Element in Engineering '\ 
Prentice-Hall, 1991, bekannt . Durch Simulation mit Hilfe Fi- 
niter Elemente werden Festigkeitsaufgaben aller Art, z. B, 
zur Spannungsverteilung oder Stabilitat, numerisch gelost. 
Beispielsweise wifd ermittelt, wie sich ein System aus star- 
ren Korpem unter auSeren Belastungen verformt und verbiegt 
und wie sich die Korper relativ zueinander verschieben, Gege 
ben ist eine rechnerverfugbare Konstruktion eines zu untersu 
chenden Systems • In der Konstruktion wird eine bestimmte Men 
ge von Punkten festgelegt, die Knotenpunkte („nodes") heiSen 
Als Finite Elemente werden die Flachen- oder Volumenelemente 
bezeichnet, die mit Hilfe der Knotenpunkte als deren Ecken 
gebildet werden, Gekrummte Flachen oder Korper, die nahe- 
rungsweise als Flachen behandelt werden, z. B. Bleche einer 
Karosserie eines Kraf tf ahrzeugs, werden hierbei oft in Scha- 
lenelemente („shell elements") zerlegt. Die Knotenpunkte bil- 
den ein Netz in der Konstruktion, weswegen der Vorgang, Kno- 
tenpunkte festzulegen und Finite Elemente zu erzeugen, Ver- 
netzung („meshing") der Konstruktion genannt wird. Je nach 
Aufgabenstellung werden die Verschiebungen dieser Knotenpunk- 
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te und/oder Rotationen der Finiten Elemente in diesen Knoten- 
punkten oder die Spannungen in diesen Finiten Elementen als 
Unbekannte eingefuhrt. Gleichungen werden aufgestellt, welche 
die Verschiebungen, Rotationen oder Spannungen innerhalb ei- 
nes Finiten Elements naherungsweise beschreiben. . Weitere 
Gleichungen result ieren aus Abhangigkeiten zwischen verschie- 
denen Finiten Elementen, z. B, daraus, daS das Prinzip der 
virtuellen Arbeit in den Knotenpunkten erfullt sein mu6 und. 
die berechneten Verschiebungen stetig sein mussen und die 
Randbedingung erfiillen mussen, daS in der Realitat Klaffungen 
oder Durchdringungen nicht auftreten. 

In vielen Fallen sind derartige Gleichungen linear in den Un- 
bekannten. Die Methode der Finiten Elemente laSt sich aber e- 
benfalls im Falle nichtlinearer Gleichungen anwenden, z. B. 
fiir Gleichungen in Form von Polynomen. Insgesamt wird ein oft 
sehr umf angreiches Gleichungssystem mit den Knotenpunkt- 
Verschiebungen/ Knotenpunkt -Rotationen, Element -Spannxingen o- 
der weitere GroSen als Unbekannte aufgestellt und numerisch 
gelost. Die Losung beschreibt beispielsweise den Verformungs- 
zustand des Systems unter vorgegebenen Belastungen, Aus die- 
ser Losung lassen sich z, B. Spannungsverteilungen, Schwin- 
gungsverhalten, Beulverhalten oder Vorhersage der Lebensdauer 
ableiten. Sind z. B. die Verschiebungen und Rotationen aller 
Knotenpunkte eines Finiten Elements bestimmt, so lafit sich 
die Spannung im Element herleiten. 

Verschiedene Korper einer Konstruktion eines Systems werden 
oft unabhangig voneinander vernetzt. Beispielsweise ist das 
System Teil der Karosserie eines zu konstruierenden Kraft - 
fahrzeuges, und die Korper sind Teilsysteme, die von ver- 
schiedenen Lieferanten zeitlich parallel konstruiert werden, 
ohne dais die Vernetzungen aneinander angepalSt werden. Weil 
die K6rper unabhangig voneinander vernetzt sind, liegen die 
Knotenpunkte auf aneinander angrenzenden Oberfiachen der Kor- 
per oft nicht aufeinander, sondern sind z. B. gegeneinander 
verschoben oder geh6ren zu Finiten Elementen unterschiedli- 
Cher Grofien und unterschiedlicher Orientierungen im Raum. 
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Derartige Vernetzungen von aneinander angrenzenden K6rpem 
werden als inkotnpatible Vernetzungen bezeichnet. 

Eine realitatsnahe Finite-Elemente-Simulation muS die Wech- 
selwirkungen iind Abhangigkeiten zwischen verschiedenen Kor- 
pern, die aufgrund der aneinander angrenzenden Oberflachen 
hervorgerufen werden, berilcksichtigen. Im Falle eines Systems 
mit einem Korper und einer Schicht sind die Wechselwirkungen 
und Abhangigkeiten zwischen Korper und Schicht zu beriLicksich- 
tigen. Gewunscht werden Finite-Elemente-Simulationen, die 
diese Wechselwirkungen und Abhangigkeiten auch bei unabhangi- 
gen und daher in der Kegel inkompatiblen Vernetzungen von 
Korper und Schicht beriicksichtigen. Denn wenn eine kompatible 
Vernetzung fur die Aufstellung des Gleichungssystems und 
Durchfuhrung der Simulationen notwendig sein wiirde, konnten 
die Korper nicht unabhangig voneinander vernetzt werden. 

Ein Verfahren zur Finite-Elemente-Simulation einer Klebever- 
bindung ist aus G. Tokar: „P\mktschweiSkleber - Eigenschaf ten 
und Berechnungsmethode fur lineare Karosseriesteif igkeiten'\ 
VDI-Berichte Nr. 1559, S. 549 - 575, 2000, bekannt , Das Sys- 
tem umfafit zwei Korper und eine Schicht, welche die beiden 
Korper verbindet. Die beiden Korper sind zwei Bleche, die 
durch eine Schicht in Form einer Klebenaht verbunden werden, 
Aufgrund auSerer Belastungen treten Verschiebungen zwischen 
den Blechen und Verformungen der beiden Bleche auf . 

Ausgehend von einer rechnerverf dgbaren Konstruktion des Sys- 
tems werden die Bleche naherungsweise durch Flachen in der 
Blechmitte reprasentiert . Fur jedes Blech wird eine Menge von 
Schalenelementen in der Blechmittenebene erzeugt . Die Ver- 
schiebungen der Knotenpunkte dieser Schalenelemente fungieren 
als Unbekannte einer Finite-Elemente-Simulation. Die Vernet- 
zungen der beiden Bleche konnen unabhangig voneinander er- 
zeugt worden sein, eine kompatible Vernetzung wird nicht vor- 
ausgesetzt. Eine sogenannte Interpolationsf lache auf derjeni- 
gen Oberf lache jedes Blechs, die der Klebenaht zugewandt ist, 
wird erzeugt. Die Klebenaht wird in der Konstruktion durch 
diese Interpolationsf lachen begrenzt und ebenfalls in Finite 
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Elements mit Knotenpunkten zerlegt. Urn Wechselwirkungen zwi- 
schen einer Interpolationsf lache als Begrenzxangsf ISche der 
Klebenaht und dem angrenzenden Blech zu beriicksichtigen, wer- 
den Knotenpunkte der Schalenelemente auf die Interpolations- 
f lache abgebildet und Interpolationen in der Interpolations - 
f lache durchgef uhrt . 

Diese in G. Tokar, a.a.O./ beschriebene kontinuums-Tnechanisch 
korrekte Abbildung der Vorgange erfordert die Einfuhrung vie- 
ler Knotenpunkte und fxihrt damit zu umf angreichen Gleichungs- 
systemen mit vielen Unbekannten. Die Aufstellung und numeri- 
sche Losung dieses Gleichungssystems ist zeitaufwendig, ins- 
besondere dann, wenn das Aufstellen und das Losen wiederholt 
durchzufuhren sind. Eine ebenfalls offenbarte Vereinf achung 
ist nicht kontinuums-mechanisch, d. h, behandelt die Klebe- 
naht nicht als Kontinuum, sondern approximiert sie durch ein- 
zelne Klebepunkte. Diese Vereinf achung fuhrt in vielen Situa- 
tionen, z. B. bei Spannungsberechnungen^ zu Simulationsergeb- 
nissen, die nicht realitatsnah genug sind, 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren be- 
reitzustellen, durch das automatisch ein Gleichungssystem fur 
eine realitatsnahe Finite-Elemente-Simulation dergestalt er- 
zeugt wird, daS die Aufstellung und das numerische Losen des 
Gleichungssystems mit einem geringeren Rechenaufwand als bei 
Anwendung bekannter Verfahren verbunden sind. Die Simulation 
soil sich auf ein System, das einen Korper und eine angren- 
zende Schicht umfaSt, beziehen. 

Die Aufgabe wird durch ein Verfahren nach Anspruch 1 gelost. 
Vorteilhafte Ausgestaltungen sind in den Unteranspruchen an- 
gegeben . 

Das automatisch erzeugte Gleichungssystem bezieht sich auf 
Knotenpunkte einer Konstruktion, die einen Korper und eine 
angrenzende Schicht umfaSt. Finite Elemente filir den Korper 
werden erzeugt . Die Schicht wird automatisch in Volumenele- 
mente zerlegt, wobei geometrischen Inf ormationen uber die 
Schicht und Vorgaben fiir die Vernetzung der Schicht verwendet 
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warden. Fur mindestens einen Knotenpimkt , der zu einem Volu- 
menelement der Schicht gehort, werden ein dem Knotenpunkt 
nachstliegendes Finites Element des Korpers und ein nachst- 
liegender Punkt auf diesem Finiten Element ermittelt. Bine 
Funktion fur einen physikalischen Zusammenhang zwischen dem 
Wert, den die physikalische Grofie im nachstliegenden Punkt 
annimmt, und den Werten, den diese GroSe in den Knotenpunkten 
des nachstliegenden Finiten Elements annimmt, wird erzeugt , 

Die physikalische GroiSe ist beispielsweise eine mechanische, 
kinematische , elektrische oder thermodynamische GroSe. Ihre 
Werte konnen Skalare oder Vektoren sein. Beispielsweise ist 
der Wert der physikalischen GroSe in einem Punkt: 

- die raumliche Verschiebung des Punktes in einem dreidimen- 
sionalen Koordinatensystem, 

- die Rotationswinkel im Punkt von demjenigen Finiten Ele- 
ment, zu dem der Punkt gehort, in einem dreidimensionalen 
Koordinatensystem und 

- die Temperatur im Punkt . 

Das physikalische Verhalten der Schicht lafit sich nur dann 
realitatsnah vorhersagen, wenn die physikalischen Abhangig- 
keiten und Wechselwirkungen zwischen dem Korper und der an- 
grenzenden Schicht berucksichtigt werden. Bei der Aufstellung 
des Gleichungssystems werden daher nicht nur die Knotenpunkte 
der Finiten Elemente des Korpers und der Schicht berucksich- 
tigt, sondern auch weitere Punkte des Korpers. Einige oder 
alle dieser weiteren Punkte konnen Knotenpunkte sein. Moglich 
ist aber auch, dafi kein weiterer Punkt ein Knotenpunkt ist. 
Erf indungsgemafi werden als weitere Punkte solche Punkte in 
Finiten Elementen des Korpers ermittelt, die Knotenpunkten 
von Finiten Elementen der Schicht nachstliegend sind. Zwi- 
schen Knotenpunkten der Schicht und nachstliegenden Punkten 
der Finiten Elemente des Kdrpers bestehen physikalische, z. 
B. mechanische, kinematische oder elektrische, Abhangigkei- 
ten. Ein Beispiel fur eirie solche Abhangigkeit ist das Prin- 
zip der virtuellen Arbeit, demzufolge die Krafte und Momente 
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zwischen den Knotenpunkten auf einer Begrenzungsf lache der 
Schicht und den nachstliegenden Punkten des angrenzenden K6r- 
pers im Gleichgewicht sind. Fiir eine realitatsnahe Simulation 
mussen diese Abhangigkeiten in das zu erzeugende Gleichungs- 
system einfliefien. 

Der Knotenpunkt, fur den das nachstliegende Finite Element 
und der nachstliegende Punkt ermittelt werden, kann sowohl 
auf einer Begrenzungsf lache als auch im Inneren der Schicht 
liegen. Die Begrenzungsf lachen konnen ebene oder gekrummte 
Flachen sein. Sie konnen zueinander parallel sein oder rela- 
tiv zueinander geneigt sein. 

Eine Moglichkeit, die Abhangigkeiten zwischen Korper und 
Schicht zu berucksichtigen, ist die, die jeweils nachstlie- 
genden Punkte der Finiten Elemente des Korpers als zusatzli- 
che Knotenpunkte mit zusatzlichen Freiheitsgraden zu verwen- 
den und dadurch die Finiten Elemente zu imterteilen. Dies 
hatte aber den Nachteil, daS das zu losende Gleichungssystem 
aufgrund der zusStzlichen Freiheitsgrade zusatzliche Unbe- 
kannte hat, z. B. die Verschiebungen dieser zusatzlichen Kno- 
tenpunkte xand Rotationen des Finiten Elements in diesen zu- 
satzlichen Knotenpunkten. 

Die Erfindung zeigt einen Weg auf, die Abhangigkeiten zwi- 
schen Korper und Schicht adaquat zu beriicksichtigen, ohne daS 
zusatzliche Freiheitsgrade auftreten und dadurch das Glei- 
chungssystem aufgrund zusatzlicher Unbekannter vergrofiert 
wird. Dank der Erfindung treten keine zusatzlichen Unbekann- 
ten auf. Durch die Berucksichtigung zusatzlicher Unbekannte 
wurde das ohnehin meist sehr umfangreiche Gleichungssystem 
noch groSer werden. Erf indungsgemaS wird eine Funktion fur 
einen physikalischen Zusammenhang zwischen dem Wert, den die 
physikalische GroSe im nachstliegenden Punkt annimmt, und den 
Werten, den diese Grofie in den Knotenpunkten des nachstlie- 
genden Finiten Elements annimmt, erzeugt. Diese Funktion 
druckt den physikalischen Zusammenhang mit zureichender Ge- 
nauigkeit aus. In alien erzeugten Gleichungen wird der Wert 
der physikalischen GroSe im nachstliegenden Punkt durch die 
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Funktion ersetzt. Die Unbekannten, die aus den Freiheitsgra- 
den des nachstliegenden Punktes resultieren, treten nicht im 
Gleichungssystem auf , sondern werden nach dem Losen des Glei- 
chungssystems mit Hilfe der Funktion fur den physikalischen 
Zusammenhang berechnet. 

Ein Beispiel fur eine solche Gleichung, in welche die Funkti- 
on eingesetzt wird, ist eine Gleichung, die aus einer oben 
beschriebenen Abhangigkeit zwischen Knotenpunkten der Schicht 
und nachstliegenden Punkten der Korper resultiert und die in 
der Simulation berucksichtigt wird. 

Durch das Einsetzen wird der Wert im nachstliegenden Punkt e- 
liminiert- Eine Unbekannte wird somit eingespart. Deshalb 
tritt der Wert im ermittelten Punkt nicht mehr als Unbekannte 
im Gleichungssystem auf, sondern hangt von den Werten der 
physikalischen GroEe in den Knotenpunkten ab, die ohnehin als 
Unbekannte des Gleichungssystems auf treten. 

Der Vorteil/ daS im Gleichungssystem Unbekannte eingespart 
werden, fallt insbesondere bei umf angreichen Systemen, z, B. 
einem Teil einer Karosserie eines Kraf tf ahrzeugs, ins Ge- 
wicht. Das Gleichungssystem umfafit dann oft Hunderttausende 
Oder gar Millionen von Unbekannten. Seine Losung erfordert 
erhebliche Rechenkapazitat \md Rechenzeit. Der Vorteil fallt 
dann noch starker ins Gewicht, wenn das Gleichungssystem wah- 
rend des Konstruierens des Systems verandert wird und daher 
mehrmals geldst und ausgewertet werden muS. Dies ist z. B. 
dann erforderlich, wenn verschiedene mogliche Konstruktionen 
eines technischen Systems verglichen werden sollen oder wenn 
verschiedene KonstruktionsstSnde durchlaufen werden und fHix 
jede Variante oder jedem Konstruktionsstand eine Finite- 
Elemente-Simulation durchgefiihrt und das physikalische Ver- 
halten der Variante oder des Konstruktionsstandes vorherge- 
sagt werden soil. 

Der Vorteil der eingesparten Unbekannten fSllt auch dann be- 
sonders stark ins Gewicht, wenn Verformungen und/oder Span- 
nungen des technischen Systems unter einer zeitlich verander- 



wo 2004/109559 



PCT/EP2004/005776 



8 

lichen aufieren Belastung vorhergesagt warden sollen. In die- 
sem Fall werden durch Diskretisierung der Zeitachse viele 
Zeitpunkte oder Zeitintervalle festgelegt. Fur jeden Knoten- 
punkt und jeden dieser Zeitpunkte wird eine Unbekannte er- 
zeugt, namlich der Wert der physikalischen GroiSe an diesem 
Knotenpunkt und zu diesem Zeitpunkt. Falls beispielsweise die 
Diskretisierung 1000 Zeitpunkte liefert, so spart das erf in- 
dungsgemaSe Verfahren bereits bei einer einzigen Anwendung N 
* 1000 Unbekannte ein, wobei N die Anzahl der Freiheitsgrade 
des gesamten Systems ist • 

Der Vorteil der eingesparten Knotenpunkte fallt dann beson- 
ders stark ins Gewicht, wenn das Gleichungssystem nichtlinea- 
re Gleichungen umfaSt und daher naherungsweise gelost wird, 
namlich iterativ durch mehrf aches Losen je eines linearen 
Gleichungssystems. Oft sind mehrere Dutzend Iterationen, in 
denen je ein lineares Gleichungssystem gelost wird, erf order- 
lich, um eine zureichend genaue Naherungslosung fur ein 
nichtlineares Gleichungssystem zu ermitteln. Dank der Erf in- 
dung werden fur jeden Iterationsschritt Unbekannte einge- 
spart - 

Das erf indungsgemaSe Verfahren lafit sich auch dann anwenden, 
wenn der Korper und die angrenzende Schicht unabhangig von- 
einander vernetzt werden und daher inkompatible Vernetzungen 
aufweisen. Daher konnen der Korper und die Schicht parallel 
konstruiert werden, z. B. von unterschiedlichen Bearbeitem, 
die sich nicht liber die Vernetzung abstimmen mussen, Weil pa- 
ralleles Konstruieren und paralleles Vernetzen ermoglicht 
wird und keine Abstimmung liber die Vernetzungen erforderlich 
ist, wird Zeit eingespart und ein simultaner Produktentwurf 
ermoglicht. Das Verfahren lafit sich fiir eine Konstruktion mit 
einem Korper und einer angrenzenden Schicht, z. B, einem 
Blech mit einer Lackschicht, anwenden. Es laSt sich auch auf 
eine Konstruktion mit zwei Korpern und einer die Korper ver- 
bindenden Schicht, z. B. eine Klebefl^che oder Klebenaht oder 
Isolierschicht, anwenden. MSglich ist, die Kdrper unabhangig 
voneinander zu vernetzen und zunachst Finite-Elemente- 
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Simulationen fiir jeden Korper unabhangig von anderen Korpern 
durchzuf ilhren . FUr die erf indungsgemaSe Finite-Elemente- 
Simulation der gesamten Konstruktion lassen sich die einmal 
erzeugten Vernetzungen der einzelnen Korper wiederverwenden . 

Das erf indungsgemaSe Verfahren erfordert keine kompatible 
Vernetzung der Korper und angrenzender bzw. verbindenden 
Schicht - Insbesondere hangt die Erzeugung der Knotenpunkte 
fur die Schicht, also die Vernetzung der Schicht, nicht von 
der Vernetzung der Korper ab und laSt sich daher gut an die 
jeweilige Auf gabenstellung, die mit Hilfe der Losung des er- 
f indungsgemafi erzeugten Gleichungs systems behandelt werden 
soil, anpassen. Beispielsweise wird je nach Aufgabenstellung 
die angrenzende bzw. verbindende Schicht in viele kleine oder 
wenige groSe Finite Elemente zerlegt . Die Dicke der Schicht 
wird beriicksichtigt - auch dann, wenn die Schicht an ver- 
schiedenen Stellen unterschiedliche Dicken aufweist. Die 
Schicht wird im Gleichungs system kontinuums-mechanisch behan- 
delt. Weiterhin lassen sich mechanische Kenngrofien der 
Schicht in Gleichungen des Gleichungssystems beriicksichtigen. 
1st die Schicht z, B. eine Klebenaht, konnen mechanische 
KenngroSen des verwendeten Klebstoffs berucksichtigt werden. 
Das mechanische Verhalten der verbindenden Schicht bei Ver- 
schiebungen des Korpers relativ zur angrenzenden Schicht oder 
bei relativen Verschiebungen von zwei Korpern, die durch die 
Schicht verbunden werden, laSt sich vorhersagen. 

Ein weiterer Vorteil des Verfahrens ist der, daS es nicht er- 
forderlich ist, ein Ersatzmodell auf zustellen. Das Auf stellen 
und Validieren eines Ersatzmodells erfordert Arbeitszeit und 
ist f ehlertrachtig . Der Vorteil, ohne Ersatzmodell auszukom- 
men, fallt insbesondere dann ins Gewicht, wenn verschiedene 
Konstruktionen eines technischen Systems mittels der Finite- 
Element e-Methode verglichen werden sollen. 

Erf indungsgemaS werden Volumenelemente fur die Schicht auto- 
matisch erzeugt . Weil die Schicht in Volumenelemente zerlegt 
wird, lalSt sich ein kontinuums-mechanisches Verhalten der 
Schicht vorhersagen, indem das Gleichungs system gelost und 



V 
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die Losung analysiert wird. Die komplette Schicht wird in das 
Gleichungssystem einbezogen. Erst durch Einbeziehung der ge- 
samten Schicht wird es erm6glicht, mindestens eine oder auch 
mehrere der folgenden Berechnungen fur die Konstruktion des 
Systems vorzunehmen : 

- geometrische Verschiebungen aufgnmd auSerer Belastungen, 
. Steifigkeit, 

. lokale Spannungsverteilung in der verbindenden Schicht und 
Position der Punkte mit maximaler Spannung, 

- momentane Festigkeit und Dauerf estigkeit im laufenden Be- 
trieb bei zeitlich verander lichen auSeren Belastungen, 

- Verhalten und AusmaS von Verf ormungen des Systems bei ei- 
nem Aufprall, 

- Schwingungsverhalten und Eigenf requenzen des Systems und 

- Beul verhalten, 

- Temperaturverteilung, 

- Verteilung des akustischen Drucks, 

- elektrisches Potential, 

- Richtung und Geschwindigkeit einer drei dimensional en Str6- 
miing, z. B. eine elektrische Stromung, die einer Fliissig- 
keit Oder eines Gases. 

Durch numerisches Losen des Gleichungssystems werden die Wer- 
te der entsprechenden physikalischen GroSe in den Knotenpunk- 
ten ermittelt, Der Wert der GroSe in dem mindestens einen 
Knotenpunkt wird ermittelt, indem die Werte in den Knoten- 
punkten in die Fxinktion fur den Zusammenhang eingesetzt wer- 
den. Die LSsung wird mittels bekannter Verfahren ausgewertet, 
um das gesuchte mechanische Verhalten vorherzusagen. 

Fur praktische Anwendungen, z. B. fur komplexere Teilsysteme 
eines Kraf tf ahrzeugs, ist es nicht in vertretbarer Zeit mog- 
lich, von Hand Volumenelemente derart zu erzeugen, daJS Simu- 
lationen auf Basis dieser Volumenelemente zu real i tat snahen 
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Ergebnissen fQhren. Dies gilt insbesondere dann, wenn die 
Vernetzung mit Vol umene lament en mehrmals durchzufiihren ist, 
z. B. mit Volumenelementen xinterschiedlicher Kantenlangen fur 
verschiedene physikalische Grofien. Die Erfindung zeigt einen 
Weg auf, wie die Vernetzung der Schicht mit Volumenelementen 
automatisch, schnell, systematisch, nachvollziehbar und in 
kurzer Zeit ausgefuhrt wird, 

Moglich ist, die Schicht senkrecht zur Schichtdicke in mehre- 
re Volumenelemente zu zerlegen. Palls z. B. die Schicht 0,8 
mm dick ist und vorgegeben ist, daJS die Schicht senkrecht zur 
Schichtdicke in zwei Volumenelemente zerlegt werden soil, so 
werden Volumenelemente erzeugt, die senkrecht zur Schichtdi- 
cke, also senkrecht zu den Begrenzungsf lachen der Schicht, 
eine Kantenlange von je 0,4 mm haben. Die Vernetzung der 
Schicht wird durch wenige und anschauliche Parameter gesteu- 
ert. Diese Parameter lassen sich so auswahlen, dafi die Ver- 
netzung fur die jeweilige Aufgabenstellung die besten Ergeb- 
nisse liefert, Vorzugsweise sind die Volumenelemente Quader, 
aber auch andere Formen von Volumenelementen sind mSglich. 

Vorzugsweise wird das erf indungsgemafie Verfahren fur jeden 
Knotenpunkt der Schicht durchgefiihrt , der auf derjenigen Be- 
grenzungsf 1 ache der Schicht liegt, die an den Korper angrenzt 
(Anspruch 2) . Damit lassen sich mechanische AbhSngigkeiten 
zwischen den Knotenpunkt en auf der angrenzenden Begrenzungs- 
flache der Schicht und den jeweils nachstliegenden Punkten 
des Korpers bei der Aufstellung des Gleichungs systems beruck- 
sichtigen, ohne hierfur zusatzliche Unbekannte zu erzeugen. 
Fur eine realitatsnahe Simulation ist gerade das mechanische 
Verhalten im Ubergangsbereich zwischen Schicht und Korper zu 
berucksichtigen . 

Fur jeden Knotenpunkt auf der angrenzenden Begrenzungsf lache 
der Schicht werden die erf indungsgemSiSen Verf ahrensschritte 
durchgef lihrt , insbesondere Ermittlung von nSchstliegendem Fi- 
nitem Element und nachstliegendem Punkt auf diesem Finiten E- 
lement, die Erzeugung der Funktion und die Ersetzung des 
Werts im ermittelten Punkt mittels der Funktion. 
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Vorzugswerden dariiber hinaus Gleichungen fur mechanische Zu- 
sammenhange zwischen Knotenpunkten der Schicht und ermittel- 
ten Piinkten der Korper erzeugt und dadurch mechanische Abhan- 
gigkeiten zwischen Schicht und Korper in vorteilhaf ter Weise 
berucksichtigt (Anspruch 3) - Fur diese Abhangigkeiten werden 
beispielsweise Finite Elemente der Form „RBE2^^ verwendet. 
Diese AbhSngigkeiten resultieren oft aus kinematischen Wech- 
selwirkungen, die sich durch eine lineare Naherung beschrei- 
ben lassen. Das erf indungsgemalSe Verfahren liefert besonders 
realitatsnahe Ergebnisse, weil durch Ermittlung der Punkte in 
den Korper die Einwirkungen, welche die Schicht auf die Kor- 
per ausubt, berucksichtigt werden, ohne daS zusatzliche Unbe- 
kannte erzeugt und verwendet werden. 

In der Ausgestaltung nach Anspruch 6 wird der Korper in der 
Simulation durch eine Flache approximiert . Im Falle eines 
dunnen Blechs als Korper ist die approximierende Flache be- 
vorzugt dessen Mittelf ISche . Als Finite Elemente der Flachen 
werden Flftchenelemente erzeugt. Die FlSche und die Flachen- 
elemente treten in der Vernetzung und Simulation an Stelle 
des Korpers und der Finiten Elemente des Korpers. 

Weil Finite Elemente in Flachen in der Regel weniger Knoten- 
punkte besitzen als Finite Elemente in Korpern, spart diese 
Ausgestaltung weitere Unbekannte ein. Damit werden Glei- 
chungssysteme mit weniger Unbekannten erzeugt, die sich 
schneller und mit weniger erforderlicher Rechenkapazitat er~ 
zeugen und numerisch losen lassen. 

Fiir die automatische Zerlegung der angrenzenden Schicht in 
Volumenelemente wird gemaS Anspruch 12 mindestens eine der 
folgenden Vorgabe verwendet - moglich ist, daS mehrere Vorga- 
ben kombiniert werden: 

- Vorgabe eines Vernetzungsverf ahrens, 

- eine untere und/oder obere Schranke fflr das Verhaltnis von 
langster zu kurzester Kante eines Volumenelements, 

- Anzahl der Volumenelemente, in welche die Schicht senk- 
recht zur Schichtdicke zerlegt werden soil, 
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- eine iintere und/oder obere Schranke f\lr die Kantenlange 
eines Volumenelements parallel zur Begrenzungsf lache der 
Schicht, 

- Vorgabe durch Projektion, 

- die geometrische Form der Volumenelemente, z. B. Quader, 
Hexaeder oder Pentagone. 

Wird die Projektion als Vorgabe festgelegt, so werden die Fi- 
niten Elemente des Korpers auf die an den Korper angrenzende 
Begrenzungsf lache der Schicht projiziert. Durch die Projekti- 
on werden die Begrenzungsf lachen der Volumenelemente vorgege- 
ben. 

Falls die Schicht sich zwischen zwei Korpern befindet und die 
Finiten Elemente des einen Korpers auf die Begrenzungsf lache 
der Schicht projiziert werden^ so ist die Vernetzung der 
Schicht in der Regel nicht kompatibel mit der Vernetzung des 
anderen K6rpers. 

Die angrenzende Sicht ist gemafi Anspruch 13 

- eine Lackschicht, 

- eine Schutzschicht oder 

- eine isolierende Schicht. 

Beispielsweise ist die Schicht eine Verzinkung, die auf einem 
Blech als dem Korper aufgetragen ist, oder eine organische 
Beschichtung, z. B. mit Polymeren, die den Korper gegen Kor- 
rosion schutzti oder eine Warme- oder K^lteisolierung oder 
eine gegen Feuchtigkeit schutzende Schicht. 

Der Begriff „ angrenzende Schicht" umfafit eine Schicht, die 
den Korper von der Umgebung des Systems isoliert, z. B. eine 
Lackschicht auf einem Blech oder eine isolierende Schicht auf 
einem Kabel. Unter dem Begriff ,,angrenzende Schicht'* wird 
auch eine Schicht verstanden, die sich zwischen zwei K6rpem 
des Systems befindet und an beide Korper angrenzt. Die 
Schicht hat beispielsweise die Funktion, die beiden Korper zu 
verbinden oder einen Mindestabstand zwischen den Korpern si- 
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cherzustellen oder Kollisionen zwischen den Korpem zu damp- 
fen Beispiele fur eine derartige Schicht sind eine Klebever- 
bindung zwischen einem Blech und einem Kunstetoff teil, eine 
ieolierende Schicht zwischen zwei Kabeln und eine Guramidich- 
tung zwischen zwei Blechen. 

Im Falle einer Konstruktion mit zwei Korpem und einer 
Schicht zwischen den Korpern wird das erf indungsgemaUe Ver- 
fahren fdr jeden Knotenpunkt der Schicht durchgefuhrt , der 
auf einer Begrenzungsf lache der Schicht liegt (Anspruch 15) . 
Hierbei wird die Schicht einmal in Volumene lament e zerlegt . 
Fur jeden Knotenpunkt der Schicht, der auf einer der beiden 
Begrenzungsf lachen der Schicht liegt und an einen der beiden 
Korper angrenzt, werden folgende Verf ahrensschritte durchge- 
f uhrt : 

die Ermittlung des jeweils nachstliegenden Finiten Ele- 
ments und des nachstliegenden Punkts, 

- die Erzeugung der Funktion fur den Zusammenhang und 

- die Eliminierung des Werts im ermittelten Punkt. 

Fur einen Knotenpunkt auf einer Begrenzungsf lache wird als 
naherliegender Korper derjenige Korper ermittelt, an den die 
Begrenzungsf lache angrenzt , Fur jeden Knotenpunkt auf einer 
Begrenzungsf lache der Schicht werden die erf indungsgemaSen 
Verf ahrensschritte durchgefuhrt , insbesondere Ermittlung von 
nachstliegendem Finitem Element und nachstliegendem Punkt auf 
diesem Finiten Element, die Erzeugung der Funktion und die 
Ersetzung des Werts im ermittelten Punkt mittels der Funkti- 
on, 

Die Ausgestaltung nach Anspruch 15 zeigt einen Weg auf, ein 
Glei Chungs system fur eine Konstruktion mit zwei Korpern und 
einer Klebeverbindung zwischen diesen zu erzeugen, ohne zu~ 
satzliche Unbekannte zu erzeugen und dadurch das Gleichungs- 
system zu vergrolSern. Die Klebeverbindung wird dabei als Kon- 
tinuum behandelt und nicht naherungsweise durch einzelne Kle- 
bepxinkte ersetzt. Fur die Aufstellung des Gleichungssystems 
brauchen lediglich Konstruktionen der beiden verbundenen Kor- 
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per, ihre Anordnung zueinander in einem System, mechanische 
KenngroiSen der Klebeverbindung und einwirkende Krafte, z. B. 
durch Einspanmong in einer Auf spannvorrichtving wahrend der 
Fertigung, vorgegeben zu sein. Die Korper brauchen noch nicht 
physikalisch produziert worden zu sein, weswegen die Vorher- 
sage zu einem friihen Zeitpunkt der Produktentwicklung vorge- 
nommen werden kann. 

Die Schicht, welche die beiden Korper verbindet, ist gemafi 
Anspruch 16 z. B. eine Klebeverbindung. Klebeverbindungen 
werden z. B. im Automobilbau zunehmend angewendet, weil eine 
SchweilSverbindung technisch nicht herstellbar ist, die zu 
verbindenden Flachen fur SchweiSer oder SchweiSautomaten 
schwer zuganglich sind oder weil die SchweiEverbindung den 
auftretenden Belastungen und Kraften nicht standhalten kann. 
Eine SchweiSverbindung ist insbesondere dann oft nicht mog- 
lich Oder unwirtschaf tlich, wenn die beiden Korper aus unter- 
schiedlichen Werkstoffen gefertigt werden, z. B. Aluminium 
und Stahl oder Aluminium und Magnesium, oder wenn mindestens 
einer der Begrenzungs flachen der Korper aus Kunststoff be- 
steht . 

Vorzugsweise wird das erf indungsgemaJS erzeugte Gleichungs sys- 
tem numerisch gelost (Anspruch 17) . Dadurch werden wenigstens 
naherungsweise die Werte ermittelt, welche die physikalische 
GroSe in den Knotenpunkten der Konstruktion annimmf. Bevor- 
zugt werden anschlieSend weitere Berechnungen vorgenommen, um 
das mechanische Verhalten des Systems zu ermitteln. Falls 
beispielsweise der Wert der physikalischen Grofie in einem 
Punkt die Verschiebung des Punktes ist, so werden durch Losen 
des Gleichungssystems die Verschiebungen aller Knotenpunkte 
ermittelt. Fur jedes Finite Element wird anschlieSend in Ab- 
hangigkeit von den Verschiebungen seiner Knotenpunkte die 
Spannung oder der Spannungstensor ermittelt, denen das Finite 
Element aufgrund der Verschiebungen ausgesetzt ist. Der Span- 
nungstensor wird ermittelt und mit einem vorgegebenen Krite- 
rium bewertet. Beispielsweise wird die maximale Spannung al- 



wo 2004/109559 



PCT/EP2004/005776 



16 

ler Finiten Elemente mit einer vorgegebenen oberen Schranke 
verglichen* 

Im folgenden wird ein Ausfuhrungsbei spiel der Erfindung an- 
hand der beiliegenden Zeichnungen naher beschrieben. Dabei 
zeigen: 

Fig. 1. ein System mit einer Klebeverbindung zwischen zwei 
FlSchen; 

Fig. 2. die Ermittlung eines nachstliegenden Finiten Elements 
und eines nachstliegenden Piinktes im Falle annahernd pa- 
ralleler FlSchen; 

Fig. 3. die Uberprufung auf Parallelitat im Beispiel der 
Fig. 2; 

Fig. 4. die Ermittlung eines nachstliegenden Finiten Elements 
und eines nachstliegenden Punktes im Falle nicht paralle- 
ler Flachen; 

Fig. 5. bilineare Interpolation in einem viereckigen Finiten 
Element mit vier Knotenpunkten in einer Ebene; 

Fig. 6. die geometrische Bedeutung der transf ormierten Koor- 
dinaten s und t; 

Fig. 7. bilineare Interpolation in einem dreieckigen Finiten 
Element . 
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Die Ausfuhrungsform bezieht sich auf eine Karosserie eines 
Kraft fahrzeugs als technischem System. Die Karosserie umfaSt 
verschiedene Bleche sowie andere Korper. Gegeben ist eine 
rechnerverfugbare Konstruktion dieses Systems. Ermittelt wer- 
den soil, wie die Karosserie sich unter einer vorgegebenen 
Belastung verformt. Urn dies zu ermitteln, werden Finite Ele- 
mente erzeugt. Die raumlichen Verschiebungen aller Knoten- 
punkte aufgrund der Belastungen sind gesucht. Diese Verschie- 
bungen bilden die Unbekannten des zu erzeugenden Gleichungs- 
systems. Jeder Knotenpunkt hat sechs Freiheitsgrade . Pro Kno- 
tenpunkt werden sechs Unbekannte erzeugt, namlich eine Unbe- 
kannte pro Freiheitsgrad; 

- drei Unbekannte fur die Verschiebung des Knotenpunkts in 

y- und z-Richtung eines dreidimensionalen Euklidischen 
Koordinatensystems 

- drei weitere Unbekannte fur die Rotation, die das Finite 
Element im Knotenpunkt aufgrund der Belastung ausfuhrt, 
wobei die Rotation durch die drei Rotationswinkel um die 
drei Koordinatenachsen ausgedruckt wird. 

In Fig. 1 wird ein Ausschnitt aus dieser Konstruktion ge- 
zeigt . Dieser Ausschnitt umf afit ein Blech und einen volumen- 
haften Korper K.l. Das Blech wird in der Finite-Elemente- 
Simulation durch eine Flache F.l in der Mitte des Blechs, al- 
so durch eine Mittelf lache, approximiert . Durch die Konstruk- 
tion ist die Dicke d_l desjenigen Blechs gegeben, das durch 
F.l approximiert wird. Der K6rper K.l umfaiSt eine zum Blech 
zeigende Begrenzungsf lache F.6. F.6 besteht aus den beiden 
Teilflachen F.6a und P. 6b. Der Abstand zwischen der Flache 
F.l und dem Kdrper K.l sowie die unterschiedlichen Orientie- 
rungen des Blechs und des Korpers K.l im Raum sind in Fig. 1 
zur Verdeutlichung stark iibertrieben dargestellt. 

Das Blech und der Korper K.l werden durch eine Klebeverbin- 
dung verbunden. Diese Klebeverbindung Kl fungiert in diesem 
Beispiel als die „angrenzende Schicht". Die zu untersuchende 
Konstruktion umf afit den Korper K.l, die Flache F.l und die 
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Klebeverbindiing Kl. Diese Klebeverbindimg Kl uberdeckt nur 
einen Teil des Blechs sowie nur einen Teil des Korpers K.l 
ixnd daher nur die Teilflache P. 6a der Begrenzungsf lache F,6 
des Korpers K.l. Die Teilflache F.6b liegt auf der an die 
Klebeverbindung Kl angrenzenden Begrenzungsf lache F.6 des 
K6rpers K.l. Der iiberdeckte Teil des Blechs wird durch eine 
Teilflache F.la der approximierenden Flache F.l reprasen- 
tiert. Die Klebeverbindung Kl, die das Blech und den Korper 
K.l verbindet; fiillt den Zwischenraum zwischen den Flachen 
F.lb und F.6c vollstandig aus. Zur besseren Darstellung ist 
sind in Fig. 1 ein Abstand zwischen der Klebeverbindung Kl 
und dem Korper K.l eingetragen. In der Konstruktion fallen 
die Flachen F.6a und F.6c zusanmien. Die Flache F.lb liegt auf 
der an die Klebeverbindung Kl angrenzenden und in Fig. 1 
nicht gezeigten Oberf lache des Blechs und damit auch auf ei- 
ner Begrenzungsf lache der Klebeverbindung Kl. F.lb ist paral- 
lel zu F.6c, weil eine Klebeverbindung mit zwei parallelen 
Begrenzungsf lachen leichter zu konstruieren ist als eine mit 
nicht parallelen Begrenzungsf lachen. Die Flachen F.la und 
F.6a sind annahernd parallel. Der Abstand zwischen F.la und 
F.lb betragt 0,5*d_l (also die halbe Dicke des Blechs). 

Die beiden Flachen F.lb und F.6a sowie die Klebeverbindung Kl 
werden unabhangig voneinander vemetzt und dabei in Finite E- 
lemente zerlegt. Die Vernetzungen konnen inkompatibel zuein- 
ander sein. Durch die Vernetzung werden die beiden Flachen 
F.lb und F.6a in Flachenelemente zerlegt, die Klebeverbin- 
dung Kl in Volumenelemente, z. B. Hexaeder. 

Die Vernetzung der Klebeverbindung Kl wird automatisch ausge- 
fiihrt. Hierbei werden folgende Inf ormationen aus der rechner- 
verfugbaren Konstruktion des Systems libernommen: 

- die r&umliche Lage der beiden Flachen F.l und F.6 und 

- die Dicke des Blechs, das durch die Flache F.l approxi- 
miert wird - in diesem Beispiel hat das Blech eine ilber 
die gesamte Ausdehnung gleichbleibende Dicke, moglich ist 
aber auch eine veranderliche Dicke, 
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In diesem Beispiel betragt die Dicke der Klebeverbindung Kl 
0,8 nun. Die Dicke iind die raumliche Ausdehnung des Uberlap- 
pungsbereichs von F.l und F.6 und damit der Klebeverbindung 
Kl warden automatisch aus diesen geometrischen Inf ormationen 
Ciber die Flachen gewonnen. Moglich ist auch, statt dessen di- 
rekt die beiden begrenzenden Teilflachen F.lb und F.6b vor- 
zugeben - 

Fur die Vernetzung der Klebeverbindung Kl werden weiterhin 
folgende vorgegebene Parameter verwendet: 

- eine untere und/oder obere Schranke fiir die Kantenlange 
eines Volumenelements der Klebeverbindung Kl parallel zur 
Begrenzungsf lache der Klebeverbindung Kl, 

- die geometrische Form der Volumenelemente, z. B. Quader, 
Hexaeder oder Pentagone, 

die Anzahl von Finiten Elementen senkrecht zur Dicke der 
Klebeverbindung Kl, 

- ein Vernetzungsverf ahren, z. B. „paving^' oder „free mes- 
hing^' , 

- die Finiten Elemente des Korpers werden auf die an den 
Korper angrenzende Begrenzungsf lache der Schicht proji- 
ziert. Durch die Projektion werden die Begrenzungsf lachen 
der Volumenelemente vorgegeben und 

- die geometrische Form der Volumenelemente/ z. B. Quader, 
Hexaeder oder Pentagone, 

Moglich ist auch, anstelle einer KantenlSnge parallel zur Be- 
grenzungsf lache die Anzahl der Volumenelemente, in welche die 
Klebeverbindung Kl parallel zu ihren Begrenzungsf lachen zer- 
legt werden soli, oder die Anzahl der Volumenelemente pro 
Langeneinheit vorzugeben. 

Weiterhin ist es moglich, die Finiten Elemente der Flache F.l 
auf die an das Blech angrenzende Begrenzungsf lache F.lb der 
Klebeverbindung Kl zu projizieren. Dadurch werden die Begren- 
zungsflachen der Volumenelemente der Klebeverbindung Kl vor- 
gegeben. Die Vernetzung der Flache F.l und damit die der Kle- 
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beverbindung Kl ist in der Regel nicht kompatibel mit der 
Vemetzung der Flache F.i. Umgekehrt ist es auch moglich, die 
Finiten Elemente der Begrenzungsf lache F.6a des Korpers K.l 
auf die an den Korper K.l angrenzenden Begrenzungsf lache F.6c 
der Klebeverbindung Kl zu projizieren, 

Vorzugsweise haben alle Volumenelemente die Form von Quadem 
Oder wenigstens von Hexaedern. Moglich sind 2. B. auch Finite 
Elemente in Form von Pentagonen. In diesem Beispiel betragt 
die Anzahl der Volumenelemente senkrecht zur Dicke der Klebe- 
verbindung Kl 2. Senkrecht zur Dicke sollen also jeweils zwei 
nebeneinanderliegende Volumenelemente erzeugt werden. Stan- 
dardmaSig haben beide Volumenelemente dieselbe Kantenlange 
senkrecht zur Dicke, so daS alle Kanten senkrecht zur Dicke 
0,8 mm : 2 « 0,4 mm lang sind. Weiterhin wird in diesem Bei- 
spiel eine Kantenlange parallel zur Begrenzungsf lache von 5 
mm in ebenen Bereichen der Klebeverbindung Kl und 4 mm in ge- 
kriimmten Bereichen vorgegeben. 

Alternativ hierzu wird nicht die Kantenlange parallel zur Be- 
grenzungsf lache vorgegeben, sondern eine untere und/oder obe- 
re Schranke fur das Verhaltnis von langster zu kurzester Kan- 
te eines Volumenelements . Beispielsweise wird ein Verhaltnis 
von 10 in gekrummten und 12,5 in ebenen Bereichen der Klebe- 
verbindung Kl vorgegeben. Wie gerade dargelegt, betragt die 
kurzeste Kantenlange 0,4 mm.\Daraus wird automatisch als Lan- 
ge der librigen Kanten eines Volumenelements 0,4mm* 12,5 =5 
mm in gekrummten und 0,4 mm * 10 = 4 mm in ebenen Bereichen 
der Klebeverbindung Kl hergeleitet. 

Fig. 2 zeigt einen Hexaeder 110.1, der zur Klebeverbindung Kl 
gehort, sowie vier Flachenelemente 100.1, 100.2, 100.3 und 
100.4 der Flache F.la und damit der Flache F.l. Beschrieben 
wird im folgenden, wie das erf indungsgemaiSe Verfahren fiir den 
Knotenpunkt 201.1 dieses Hexaeders durchgefiihrt wird. 

Nachdem die Vemetzung der Konstruktion abgeschlossen ist, 
werden die physikalischen Zusammenhange und Randbedingungen 
erganzt . Ein Beispiel fiir einen solchen Zusammenhang be- 
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schreibt die Spannung in einem Finiten Element abhSngig von 
der Verschiebung seiner Knotenpunkte . Abhangig von der Ver- 
schiebung der Knotenpunkte wird ein Dehnungstensor e des Fi- 
niten Elements bestimmt. Vorgegeben ist eine Steif igkeits- 
Matrix (/^compliance matrix^^) D. Zwischen dem Spannungstensor 
a des' Finiten Elements iind dem Dehnungstensor s besteht der 
Zusammenhang a « D * s. 

Moglich ist, daS die Verformungen aus einer Temperaturveran- 
derung AT resultieren. Sei a der Ausdehnungs-Koef f izient des 
f-Oir die Fertigung des jeweiligen Korpers verwendeten Werk- 
stoffs. Dann besteht der Zusammenhang a = D * (s - a* AT) . 

Weiterhin wird der Zusammenhang zwischen einwirkender Kraft F 
und Verformung U bestimmt* Aus Eigenschaf ten der Werkstoffe, 
die fiir die Herstellung des jeweiligen Korpers verwendet wer- 
den, ,z. B. Elastizitats-Modul und Poisson-Zahl , und aus der 
Geometrie des Korpers wird eine Steif igkeits-Matrix K des 
Korpers hergeleitet . Zwischen der Verformung und der einwir- 
kenden Kraft besteht der Zusammenhang U = K * F. Moglich ist, 
dafi einige Komponenten von U bekannt sind, z, B. gleich Null 
sein mussen, und einige Komponenten von F bekannt und andere 
unbekannt sind. 

Das zum Knotenpunkt 201.1 nachstliegende Flachenelement wird 
ermittelt. Zunachst werden die Begrenzungsf lachen des Hexa- 
eders 110.1 ermittelt, die auf eine Flache weisen, die an die 
Klebeverbindung Kl angrenzt. Im Falle des Hexaeders 110.1 
werden die Begrenziingsf lichen 111.1 und 111.2 ermittelt. Ge- 
sucht sind das zum Knotenpunkt 201.1 nachstliegende Finite E- 
lement, das die Form eines Flachenelements hat, und der 
nachstliegende Punkt auf diesem Flachenelement. 

Ein schnell durchzufuhrendes und durch Fig, 2 illustriertes 
Verfahren sieht vor, eine Normale 210.1 auf die Begrenzungs- 
flache 111.1 zu erzeugen, die durch den Knotenpunkt 201.1 
veriauft. Diese Normale 210.1 trifft die Flache F.la im Fla- 
chenelement 100.2 und dort im Punkt 205.1. Dieser Schnitt- 
punkt 205.1 wird als nachstliegender Punkt ermittelt. Das 
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Fiachenelement 100.2, zu dem der Schnittpunkt gehort, wird 
als nachstliegendes Finites Element ermittelt. 

Die Normale auf die Begrenzungsf lache 111,1 wird nacheinander 
in die Knotenpunkte des Hexaeders 110.1 verschoben. In Fig, 2 
ist zusatzlich die Normale 210.3 durch den Knotenpunkt 201.2 
gezeigt. Diese schneidet die Flache Fl.a im Punkt 205.2. Fiir 
den Knotenpunkt 201.2 werden als nachstliegender Punkt der 
Punkt 205.2 und als nachstliegendes Finites Element das Fla- 
chenelement 100.1 ermittelt. 

Falls die Begrenzungsf lache 111.1 und die Flache F.la paral- 
lel zueinander liegen, so ist der Schnittpunkt 205.1 von 
210.1 und 100.2 genau der nachstliegende Punkt. Falls der in 
Fig. 3 gezeigte Winkel 23 0.1 zwischen der Begrenzungsf lache 

111.1 und die Flache F.la klein ist, so liegt der Schnitt- 
punkt 205.1 nahe dem nachstliegenden Punkt. In einer Ausfiih- 
rungsform wird daher zusatzlich gepriift, ob der Winkel 230.1 
zureichend klein ist. Diese Prdfung illustriert Fig. 3. Zu- 
satzlich zur Normalen 210.1 durch 201.1 wird eine Normale 

210.2 durch dasjenige Flachenelement 100.2 erzeugt, die im 
Beispiel der Fig. 3 als nachstliegendes Finites Element er- 
mittelt wurde. Der Winkel 230.1 zwischen diesen beiden Norma- 
len wird ermittelt. Im Falle paralleler Flachen betragt der 
Winkel 0. Ist dieser Winkel kleiner oder gleich einer vorge- 
gebenen oberen Schranke, so werden der Schnittpunkt 2 05.1 und 
das Flachenelement 100.2, die wie oben beschrieben ermittelt 
wurden, als. nachstliegender Punkt bzw. nachstliegendes Fini- 
tes Element verwendet. 

Ist der Winkel 230.1 grofier als die obere Schranke, so wird 
bevorzugt das im folgenden beschriebene Verfahren durchge- 
fuhrt. Zunachst wird der zu 201.1 nachstliegende Knotenpunkt 
in der Flache F.la ermittelt. Im Beispiel der Fig. 4 ist dies 
der Knotenpunkt 200.1. Um die Abstandsberechnungen durchzu- 
fiihren, werden vorab alle Abstande zwischen verschiedenen 
Knotenpunkten der Konstruktion ermittelt und in einer Tabelle 
zwischengespeichert. Der Knotenpunkt gehort zu den vier FIS- 
chenelementen 100.1, 100.2, 100.3 und 100.4. Liegen diese 



wo 2004/109559 



PCT/EP2004/005776 



23 

vier Flachenelemente in einer Ebene, so wird eine Normale 
210,4 auf den vier Flachenelementen durch 201-1 ermittelt. 
Diese Normale schneidet die Flache Fl.a im Punkt 205.3, der 
zum Flachenelement 100,4 gehort . Damit sind 205,3 und 100.4 
ermittelt. Falls die vier Flachenelemente nicht in einer Ebe- 
ne liegen, werden vier Normalen und vier Schnittpunkte ermit- 
telt und der nachstliegende Schnittpunkt als nachstliegender 
Punkt verwendet. 

Fig. 5 veranschaulicht die bilineare Interpolation fur ein 
beliebiges viereckiges Flachenelement. Das Flachenelement hat 
die vier Knotenpunkte 200. a, 200. b, 200. c und 200. d, die in 
einer Ebene liegen, als Ecken. In dem Gleichungssystem, das 
mit Hilfe der Finiten Elemente erzeugt wird, treten die Ver- 
schiebungen P. a, P.b, P.c und P.d dieser vier Knotenpunkte 
als vier Unbekannte auf. Als nachstliegender Punkte wurde der 
Punkt 205'. X des Flachenelements ermittelt. Die Verschiebung 
P.x des Punktes 205.x wird naherungsweise als lineare Funkti- 
on der Verschiebungen P. a, P.b, P.c und P.d ausgedruckt : 

P.x = g.a * P-a + g.b * P.b + g.c * P.c + g.d * P.d 

In Gleichungen des zu erzeugenden Gleichungs systems wird P.x 
durch die Linearkombination auf der rechten Seite dieser li- 
nearen Funktion ersetzt. 

Die Gewichtsf aktoren g.a, g.b, g.c und g.d werden so be- 
stimmt, daS gilt: g.a + g.b + g.c + g.d = 1. Weiterhin werden 
sie so bestimmt, daS sie sich linear in den Koordinaten von 
205.x verandern, wenn 205.x im Flachenelement bewegt wird, 
\ind dais gilt : 

- Wenn 205.x = 200. a, dann ist P.x = P. a und g.a = 1 und g.b 
= g.c = g.d = 0. 

- Wenn 2 05.x = 2 00.b, dann ist P.x = P.b und g.b = 1 und g.a 
= g.c = g.d = 0. 

- Wenn 205.x = 200.C, dann ist P.x = P.c und g.c = 1 und g.b 
= g.c = g.d = 0, 
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- Wenn 205.x = 200-d, dann ist P.x = P.d und g.d = 1 und g.b 
= g-c = g-d = 0. 

In Fig. 5 sind die vier Mittelpunkte 200. ab, 200,bc, 200. cd 
und 200. ad der vier Seiten des FlSchenelements eingetragen. 
Die Verbindungsstrecken 210. e von 200. ad nach 200. be und 
210 -f von 200. ac nach 200. bd schneiden sich in 200. m. .200. m 
liegt auf der Mitte zwischen 200. ad und 200. be und ebenfalls 
auf der Mitte zwischen 200. ab und 200. cd. Die beiden Verbin- 
dungsstrecken 210. e und 210. f halbieren sich also. 

Die Lage des Pxinktes 205, x im Flacheneleraent vor der Ver- 
schiebung wird durch zwei Koordinaten s und t beztiglich eines 
Koordinatensystetns, dessen Achsen die beiden Verbindungsstre- 
cken 210. e und 210. f sind, beschrieben. Diese Koordinaten s 
und t werden so festgelegt, daS folgendes erfullt ist: 

- Falls 205.x - wie in Fig. 5 - zwischen den Punkten 200. a, 
200. ab, 200. ad und 200. m liegt, ist -1 <= s <= 0 und -1 <= 
t <= 0. 

- Falls 205.x zwischen den Punkten 200. ab, 200. b, 200. be und 
200. m liegt, ist 0 <= s 1 und -1 <= t <« 0. 

- Falls 205.x zwischen den Punkten 200. be, 200. c, 200. cd und 
200. m liegt, ist 0 <= s <= 1 und 0 <« t <= 1. 

- Falls 205.x zwischen den Pxinkten 200. cd, 200. d,. 200. ad und 
200. m liegt, ist -1 <= s <= 0 und 0 <= t <^ 1. 

- Fur alle Punkte auf der Strecke von 200. d nach 200. a gilt: 
s = -1 und -1 <= t <= +1. 

- Fur alle Punkte auf der Strecke von 2 00. a nach 200.b gilt: 
t = -1 und -1 <= s <« +1. 

- Fur alle Punkte auf der Strecke von 200. b nach 200. c gilt: 
s = 1 und -1 <« t <= +1. 

- Fur alle Punkte auf der Strecke von 200.c nach 20O.d gilt: 
t = 1 und -1 <= s <= +1. 

- Fixr alle Punkte auf der Strecke von 2 0O.cd nach 200.ab 
gilt: s = 0 und -1 <= t +1. 
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- Flir alle Punkte auf der Strecke von 200. ad nach 200. be 
gilt: t = 0 und -1 <= s <= +1. 

- Flir den Knotenpunkt 200. a gilt: s = -1 und t = -1. 

- Fur den Knotenpunkt 200. b gilt: s = +1 und t = -1. 

- Fur den Knotenpunkt 200. c gilt: s = +1 und t = +1. 

- Fur den Knotenpunkt 200. d gilt: s = -1 und t = +1. 

Die vier Gewichtsf aktoren werden durch folgende Rechenvor- 
schrift bestimmt: 

g.a = M * (1-s) * (1-t) 

g.b = M * (1+b) * (l-t) 

g.c = M * (1+s) * (1+t) 

g.d = ^ * (1-s) * (1+t) 

Fig. 6 illustriert die geometrische Bedeutung der beiden Ko- 
ordinaten s und t. Die Lage jedes Punkts des in Fig. 5 ge- 
zeigte Vierecks wird nach einer Koordinatentransf ormation 
durch zwei Koordinaten (x, y) beschrieben. In Fig. 6 im lin- 
ken Viereck ist der Punkt 205.x durch die Koordinaten (x, y) 
beschrieben, der Schnittpunkt 200. m durch (x_ni, y_rn) , die E- 
cken 200. a, 200, b, 200. c und 200. d durch die Koordinaten 
(x_a, y_a) , (x_b, y_b) , (x_c, y_c) , (x_d, y_d) . Das Quadrat 
auf der rechten Halfte von Fig. 6 veranschaulicht die Berech- 
nung der Koordinaten s und t. Der Punkt (x, y) hat nach der 
Transformation in das natiirliche Koordinatensystem im Quadrat 
die Koordinaten (s, t) . Der Punkt (x_m, y_m) hat die Koordi- 
naten (0, 0), der Punkt (x_a, y_a) (-1, -1), der Punkt (x_b, 
y_b) (1, -1), der Punkt (x_c, y_c) (-1, 1) und der Punkt 
(x_d, y_d) (-1. 1) . 

Die Koordinaten s und t werden bevorzugt gemaS der folgenden 
Berechnungsvorschrift berechnet: 

Falls der Punkt (x, y) auf der Kante zwischen (x_a, y_a) und 
(x_b, y_b) liegt, ist 
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2x-x a-x_b _ 

s = = ==- una 

xj>-x_a 

t - -1. 

Falls der Punkt (x, y) auf der Kante zwischen (x_c, y_c) und 
(x_d, y_d) liegt, so gilt 

2x-x c-x d _ 

s = == — vmd 

x_c-x_d 

Falls der Punkt (x, y) auf der Kante zwischen (x_a, y_a) und 
(x_d, y_d) liegt, so gilt 

s = -1 und 
t = -1. 

Falls der Punkt (x, y) auf der Kante zwischen (x_b, y_b) und 
(x__c, y_c) liegt, so gilt 

s = 1 und 
2y-yJ>-y_c 



t = - 



y_c-y_b 



Falls der Punkt (x, y) im Inneren des Vierecks von Fig. 6 
links liegt, ist 

^_-Bi_-^^ Oder t= Bi + -n/A ^ nachdem, welche Berechnungs- 

2Ai 2Ai 
vorschrift zu einem Wert von t zwischen 0 und 1 fiihrt, und 

4x-Ax-Bxt 
s = — -z — — 
Cx+Dxt 

Hierbei sind 
A^B)^ -4A\C\ 

Ai = ByDx - DyBx 

Bl = -4yDx + AyDx + ByCx - CyBx + 4xDy " AxDy 
Cl = -4yCx + AyCx + 4xCy " CyAx 
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Ax 




x_a + 


x_b + x_c + 


x_d 


Bx 




-x_a 


- x_b + x_c 


+ x_d 


Cx 




-x_a 


+ x_b + x_c 




Dx 




x_a - 


x__b + x_c - 


x_d 


Ay 




y_a + 


y_b + y_c + 




By 




-y_a 


- y_b + Yjc 


+ y_d 


Cy 




-y_a 


+ y_t> + 




By 




y_a - 


y_b + y_c - 


y^d 



Sobald also die Koordinaten s und t und die vier Verschiebun- 
gen P. a, P.b, P-c und P.d bestitnmt sind, wird die Verschie- 
bung P.x durch Einsetzen der dann bekannten Koordinaten und 
Verschiebungen bestimmt • 

Analog wird der Fall eines Flachenelements mit drei Knoten- 
punkten 200. a, 200. b und 200. c behandelt. Diesen Fall illust- 
riert Fig. 7. Die Verschiebung P.x des Punkts 205.x wird na- 
herungsweise durch die Funktion 

P.x = g.a * P. a + g.b * P.b + g.c * P-c 

ausgedriickt. Die Gewichtsf aktoren g.a, g.b und g.c werden so 
bestimmt, daS gilt: g.a + g.b + g.c = 1. Die Lage des Punktes 
205.x im Flachenelemente bezuglich eines- Koordinatensystems , 
dessen Achsen die Strecke von 200. a nach 200. b und die Stre- 
cke von 200. a nach 200. c sind, wird durch zwei Koordinaten s 
\ind t festgelegt. Die Koordinaten s und t werden so festge- 
legt, dafi folgendes gilt: 

0 <= s <« 1 und 0 <= t <= 1 und s + t <= 1. 

Diese Koordinaten s und t werden so festgelegt, daS folgendes 
erfullt ist: 

. Falls 205.x auf der Strecke von 200. a nach 200. b liegt, 
ist s+t = 1. 

- Falls 205.x auf der Strecke von 200. a nach 200. c liegt, 
ist t = 0. 
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- Falls 205.x auf der Strecke von 200 -b nach 200. c liegt, 
ist s = 0. 

" Falls 205.x « 200. a ist, ist s = 1 und t = 0. 

- Falls 205 -X - 200. b ist, ist s = 0 und t = 1. 

- Falls 205. X == 200. c ist, ist s = 0 und t = 0. 

Die Gewichtungen g.a, g.b und g.c werden abhangig von s und t 
wie folgt festgelegt: 

g.a = s, 

g.b = t und 

g.c=l-s-t. 

Das Dreieck wird wie oben beschrieben durch eine geeignete 
Koordinatentrans format ion so transf ormiert , dafi der Punkt 
205.x die Koordinaten (x, y) hat. Der Punkt 20O.a hat die Ko- 
ordinaten (x_a, y_a) , der Punkt 200. b die Koordinaten (x_b, 
y_b) und 2 00.c die Koordinaten (x_c, y_c) . Die Koordinaten s 
und t werden wie folgt bestimmt: 

Falls (x, y) auf der Kante zwischen (x_a, y_a) und (x_c, y_c) 
liegt, ist 

s = ^ und 

y_a-y_c 

t = 0 

Falls (x, y) auf der Kante zwischen (x_b, y_b) und (x_c, y_c) 
liegt, ist s = 0 xand 

y_b-y_c 

Falls (x, y) auf der Kante zwischen (x_a, y_a) und (x_b, y_b) 
Oder im Inneren des Dreiecks liegt, ist 




_t(x b-x c) + (x c-x) 
x_a-x_c 

Hierbei ist 
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(y a-y c)(x b-x c) . , . 
A = — ^" ~ — —-^ (yj) - y^c) 
x_a-x_c 

und 

^ (y a-y c)(x c-x) . . 

B = — " (y_c - y) 

x_a-x_c 

Falls ein Flachene lament vier Ecken hat, die nicht in einer 
Ebene liegen, so wird einer der folgenden beiden Schritte 
ausgef uhrt : 

Die Koordinaten s und t werden nicht bezuglich des Flachen- 
elements mit den Ecken 200. a, 200. b, 200. c und 200. d be- 
stimmt, sondern bezuglich eines approximierenden Flachenele- 
ments mit den vier in einer Ebene liegenden Ecken 200. a', 
200, b\ 200-c' und 200. d^. Die Ecken werden z. B. so be- 
stimmt, daS 200. a^ = 200. a, 200. b' = 200. b und 200. c^ = 200. c 
gilt. 200. d' ist der FujSpunkt einer Normal en durch 200. d auf 
der Ebene, die durch 200. a, 200. b und 200. c bestimmt ist. O- 
der die vier Ecken werden so bestimmt, dafi 

II 200. a^ ' 200. a || + 1| 200. b^ - 200. b || + 
II 200. c^ - 200. c II + II 200. d^ - 200. d || 

minimiert wird. Hierbei ist 200. a' - 200. a der Vektor von 
200. a' nach 200. a und || x || die Euklidische Lange eines 
Vektors x. 

Der' zweite mogliche Schritt sieht vor, daJS automatisch ent- 
schieden wird, ob 205.x im Dreieck mit den Ecken 200. a, 200. b 
und 200. c Oder im Dreieck mit den Ecken 2 00. a, 200. b und 
200.d liegt. Fiir das Dreieck, in dem 205.x liegt, wird das o- 
ben beschriebene Verfahren fiir dreieckige FlSchenelemente 
durchgefOhrt . 

Software -Programme zur Vernetzung von Konstruktionen besitzen 
eine elektronische Bibliothek mit Bausteinen. In dieser Bib- 
liothek ist typischerweise ein Baustein namens RBB3 (RBE = 
Rigid Body Element'') enthalten, das die Anwendung der Ge- 
wichtsfaktoren g.a, g.b, g.c und g.d erm6glicht. 
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Nachdem die Vernetzung der Begrenzungsf lache F.6 des Korpers 
K.l, der Mittelf lache F.l des Blechs und der verbindenden der 
Klebeverbindung Kl abgeschlossen sind und das Gleichungs sys- 
tem erzeugt worden ist, wird das Gleichvmgssystetn mit einem 
kommerziellen Sof tware-Werkzeug fur die Finite-Elemente- 
Methode (FEM-Werkzeug) geldst. 

Der Fachmann kennt verschiedene FEM-Werkzeuge, z. B. 

- MSC.NASTRAN und MSC.PATRAN, beide beschrieben unter 
http : //www.mscsof tware . com/products/ , abgefragt am 5. 2. 
2003, 

- ABAQUS, beschrieben unter 

http ; / /www , hks . com/product s /product s overview . html , abge- 
fragt am 5, 2. 2 003, 

PAMCPASH fur Finite -Element e-Simulationen von Kollisionen, 
beschrieben unter http://www.esi- 

group-com/products/crash/index,php, abgefragt am 5. 2. 
2003. 

Die Losung liefert fur jeden Knotenpunkt der Konstruktion den 
Wert, den die physikalische Grofie in diesem Knotenpunkt an- 
nimmt. Durch Einsetzen in die Funktion werden die Werte der 
physikalischen Grofie in den ermittelten nachstliegenden Punk- 
ten berechnet. Die Losung wird ausgewertet, urn die Konstruk- 
tion des Systems zu analysieren- 
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1. Verfahren zum automatischen Aufstellen eines Gleichungs- 
systems zur Beschreibung eines physikalischen Verhaltens 
eines technischen Systems, 

wobei eine rechnerverf ugbare Konstruktion vorgegeben ist, 
die einen Korper und eine Schicht (Kl) , die an den Korper 
angr en z t , umf aSt , 

mit den Schritten 

- Erzeugen von Finiten Elementen (100.1, 100.2, ...) fur 
den Korper, 

- vollstandiges Zerlegen der Schicht (Kl) in Finite Ele- 
mente (110.1), die Volumenelemente sind und automa- 
tisch unter Verwendung von geometrischen Inf ormationen 
uber die Schicht (Kl) und Vorgaben fur die Vernetzung 
der Schicht (Kl) erzeugt werden, und 

- Aufstellen eines Gleichungssystems gemafi der Finite- 
Elemente-Methode, in dem als Unbekannte die Werte auf- 
treten, die eine raumlich veranderliche physikalische 
GroSe in den Knotenpunkten (200.1, 200.2, 201.1, 
201.2, ...) der Finiten Elemente des Korpers und der 
Schicht (Kl) annimmt, 

- wobei fur mindestens einen Knotenpunkt (201.1, 201.2), 
der zu einem Volumenelement (110.1) der Schicht (Kl) 
gehort , 
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- ein dem Knotenpunkt nachstliegendes Finites Element 
(100.1, 100.2, ...) des Korpers und ein nSchstlie- 
gender Punkt (205.1, 205.2, ...) auf diesem Finiten 
Element ermittelt werden, 

- eine Fimktion f\ir einen physikalischen Zusammenhang 
zwischen dem Wert, den die physikalische GroSe im 
nachstliegenden Punkt (205.1, 205.2, ...) annimmt, 
und den Werten, den diese GroSe in den Knotenpunkten 
(200.1, 200.2, ...) des nachstliegenden Finiten Ele- 
ments (100.1, 100.2, ...) annimmt, erzeugt wird 

• und beim Aufstellen des Gleichungssystems der Wert 
der physikalischen GroSe im nachstliegenden Punkt 
(205.1, 205.2, ...) durch Einsetzen der Funktion e- 
liminiert wird. 

2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, 

daS fur jeden Knotenpunkt (201.1, 201,2) der Schicht 
(Kl) , der auf der an den Korper angrenzenden Begrenzungs- 
flache (Fl.b) der Schicht (Kl) liegt, 

- die ErmittlxHig des jeweils nachstliegenden Finiten E- 
lements (100.1, 100.2, ...) und des nachstliegenden 
Punk'ts (205.1, 20'5.2, .'..)",• 

- die Erzeugung der Funktion fiir den Zusammenhang 

- und die Eliminierung des Werts im nachstliegenden 
Punkt (205.1, 205.2, ) 

durchgefiihrt werden. 

3. Verfahren nach Anspruch 2, 

dadurch gekennzeichnet, 

dafi Gleichungen fur physikalische Zusammenhange zwischen 
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- den Werten, welche die physikalische GroSe in einer 
Menge von Knotenpunkten (201.1, 201.2), die auf der 
angrenzenden Begrenzungsf lache (Fl.b) det Schicht (Kl) 
1 iegen , annimmt 

- und den Werten, welche die physikalische GroSe in ei- 
ner Menge von Punkten des Korpers, die jeweils zu ei- 
nem Knotenpunkt (201.1, 201.2) der Menge nSchstliegend 
sind, annimmt, 

erzeugt und beim Aufstellen des Gleichungssystems verwen- - 
det werden. 

4. Verfahren nach einem der Anspruche 1 bis 3, 
dadurch gekennzeichnet, 
dais 

- mindestens ein Finites Element (110. l) der Schicht 
(Kl) erzeugt wird, das einen Knotenpunkt im Inneren 
der Schicht (Kl) besitzt, 

- und eine Gleichung fur den Zusammenhang zwischen den 
Werten, die eine physikalische GroEe in den Knoten- 
punkten (201.1, 201-2) dieses Finiten Elements (110.1) 
der Schicht (Kl) annimmt, erzeugt und beim Aufstellen 
des Gleichungssystems verwendet wird. 

5. Verfahren nach einem der Anspruche 1 bis 4, 
dadurch gekennzeichnet, 

dag der Wert der physikalischen Gr66e in einem Punkt min- 
destens einen der folgenden Bestandteile umf afit : 

- die raumliche Verschiebung des Punktes in einem drei- 
dimensionalen Koordinatensystem, 

- die Rotationswinkel im Punkt von demjenigen Finiten E- 
lement, zu dem der Punkt gehort, in einem dreidimensi- 
onalen Koordinatensystem, 
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- die Temperatur im Punkt, 

- die Temperaturveranderung im Punkt, 

- das elektrische Potential im Punkt, 

- der akustische Druck im Punkt, 

- der Richtungsvektor, den eine Stromung im Punkt hat, 

- die Geschwindigkeit, den eine Str6mung im Punkt hat. 

6. Verfahren nach einem der Anspruche 1 bis 5, 
dadurch gekennzeichnet, 
dafi 

- eine Flache (F.l), die zum Korper gehort, erzeugt 
wird, 

- Finite Elemente (100.1, 100.2, ...) fiir die Flache 
(F.l) erzeugt werden 

- und als nachstliegendes Finites Element ein nachstlie- 
gendes Finites Element (100.1, 100.2, ...) der Flache 
(F.l) ermittelt wird. 

7. Verfahren nach Anspruch 6, 

dadurch gekennzeichnet, 

daS die Funktion fiir den Zusammenhang durch bilineare In- 
terpolation uber den Knotenpunkten (200.1, 200.2, ,.,) 
des nachstliegenden Finiten Elements (100.1, 100.2, ...) 
erzeugt wird. 

8 . Verfahren nach Anspruch 6 oder Anspruch 7 , 
dadurch gekennzeichnet, 
daS 

- eine Normal e (210.1, 210.3) durch den Knotenpunkt der 
Schicht (Kl) auf einer dem Korper nachstliegenden Be- 
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grenzungsflache (111.1) des Volumenelements (110.1) 
erzeugt wird, 

- der Schnittpunkt (205.1, 205,2, ...) dieser Normalen 
(210.1, 210.3) mit der Flache (P.l) ermittelt und als 
nachstliegender Punkt verwendet wird 

- und ein Finites Element (100.1, 100.2, ...) der Flache 
(F.la), zu dem der Schnittpunkt (205.1, 205.2, ...) 
geh6rt, als nachstliegendes Finites Element verwendet 
wird. 

9. Verfahren nach einem der Anspruche 6 bis 8, 
dadurch gekennzeichnet, 
daS 

- eine Normale (210.4) durch den Knotenpunkt der Schicht 
(Kl) auf das nachstliegende Finite Element (100.4) er- 
zeugt wird 

- und als nachstliegender Punkt der FuiSpunkt (205.3) der 
Normalen (210.4) im Finiten Element (100.4) ermittelt 
wird. 

10. Verfahren nach einem der Ansprdche 6 bis 9, 
dadur'ch gekennzeichnet, 
dafi der Korper ein Blech ist 

und die Flache (F.l) in der Mitte des Blechs erzeugt 
wird. 

11. Verfahren nach einem der Anspriiche 6 bis 10, 
dadurch gekennzeichnet, 
dafi 

- die Finiten Elemente (100.1, 100.2, .,.) der Flache 
(F. 1) Schalenelemente 
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- und die Volumeneletnente (110.1) der Schicht (Kl) Qua- 
der 

sind. 

12. Verfahren nach einem der Anspriiche 1 bis 11, 

dadurch gekennzeichnet, 

daS ffir die Erzeug\ing der Volumenelemente (110.1) in der 
Schicht (Kl) mindestens eine der folgenden Vorgaben ver- 
wendet werden: 

Vorgabe eiries Vernetzungsverf ahrens, 

- eine untere und/oder obere Schranke fur das Verhaltnis 
von langster zu kiirzester Kante eines Volumenelements, 

- Anzahl der Volumenelemente, in welche die Schicht 
senkrecht zur Schichtdicke zerlegt werden soil, 

. eine untere und/oder obere Schranke fUr die Kantenlan- 
ge eines Volumenelements parallel zur Begrenzungsf la- 
che der Schicht, 

Projektion der Finiten Elemente des Korpers auf die an 
den K6rper angrenzende Begrenzungsf lache der Schicht 
und 

" die geometrische Form der Volumenelemente. 

13. Verfahren nach einem der Anspruche 1 bis 12, 
dadurch gekennzeichnet, 
daS die angrenzende Schicht (Kl) 

- eine Lackschicht, 

eine Schutzschicht oder 
eine isolierende Schicht 
auf dem Kdrper ist , 
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14. Verfahren nach einem der Anspruche 1 bis 13, 
dadurch gekennzeichnet, 
daS 

- das System einen weiteren Korper (K,l) umfaSt, 

- die Schicht (Kl) an beide Korper angrenzt, 

- finite Elemente fxir den weiteren Korper (K.l) erzeugt 
werden 

. der dem Knotenpunkt naherliegende Korper des Systems 
ermittelt wird 

- iind als nachstliegendes Finites Element ein nachstlie*- 
gendes Finites Element des naherliegenden Korpers er- 
mittelt wird- 

15. Verfahren nach Anspruch 14, 

dadurch gekennzeichnet, 

daS fiir jeden Knotenpunkt der Schicht (Kl) , der auf einer 
der beiden Begrenzungsf lachen der Schicht (Kl) liegt und 
an einen der beiden K6rper angrenzt, 

- die Ermittlung des jeweils nachstliegenden Finiten E~ 
lements und des nachstliegenden Punkts, 

- die Erzeugxmg der Funktion fur den Zusammenhang 

- und die Eliminierung des Werts im ermittelt en Punkt 
durchgefilhrt werden. 

16. Verfahren nach Anspruch 14 oder Anspruch 15, 
dadurch gekennzeichnet, 
daS die Schicht (Kl) 

« eine Klebeverbindung, 

- eine Abdichtung oder 
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- eine Isolierung 

zwischen den beiden Korpern des Systems ist. 

17. Verfahren nach einem der Anspriiche 1 bis 12, 

dadurch gekennzeichnet, 

daJS die Werte, welche die physikalische GrolSe in den Kno- 
tenpunkten der Konstruktion annimmt, ermittelt werden, 

indem das Gleichungssystem numerisch gelost wird. 

18. Computerprogramm-Prodiikt, das direkt in den internen 
Speicher eines Computers geladen werden kann imd Soft- 
wareabschnitte umfaSt, mit denen ein Verfahren nach einem 
der Anspriiche 1 bis 17 ausgefiihrt werden kann, wenn das 
Produkt auf einem Computer lauft. 

19. Computerprogramm- Produkt, das auf einem von einem Compu- 
ter lesbaren Medium gespeichert ist und das von einem 
Computer lesbare Programm-Mittel aufweist, die den Compu- 
ter veranlassen, ein Verfahren nach einem der AnsprQche 1 
bis 17 auszufuhren. 
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